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ABSTRACT 

We present the rest- frame optical spectrum of a strongly lensed galaxy at redshift z = 1.7 behind the clus- 
ter Abell 1689. We detect the temperature sensitive auroral line [Olll] A4363, which allows the first direct 
metallicity measurement for galaxies at z > 1. Our high signal-to-noise spectrum indicates that the target is an 
extremely low metallicity star-forming galaxy, with 12 + log(0/H) = T.S^q ^ from the direct Tg-based method. 
We estimate an intrinsic absolute B-hand magnitude of Mb = — 18.3±0.1, with a stellar mass of 4.4± 1.2x10^ 
M0 . This galaxy extends the luminosity — metallicity relation of star-forming galaxies at z > 2 by more than 
an order of magnitude. Given the double-nuclei like morphology and the velocity profile of Ha, we tentatively 
suggest that it could be a merger or a proto-rotating disk galaxy. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift 
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1. INTRODUCTION 

Five years ago, our knowledge about galaxy evolution still 
had a glaring gap between redshifts 1.5 and 3, termed the 
redshift desert. The desert is even more desolate for galaxy 
chemical evolution studies, where high signal-to-noise (S/N) 
spectra are needed for reliable measurements. However, great 
progress has been made in the last few years in studying 
galaxies at z~ 2 and we are beginning to harvest the "desert". 
To date, metallicities are available for a few tens of individual 
galax ies at 1.5 < z < 3 (.Shaolev et al1l2004l: lErb et al.l l2003[ 
l2QQ6l:lHavashi et al.ll2QQ9ir These galaxies are selected using 
broad-band c olors either in the UV (Lyman Break technique; 
ISteidel etal.l 11996. 200 1) or using B, z, and i^-band colors 
(BzK selection; Daddi et al.l l2004|) . The Lyman break and 
BzK selection techniques favor the galaxies that are lumi- 
nous in the UV or blue and could therefore be biased against 
low luminosity (low-metallicity) galaxies, and dusty (poten- 
tially metal-rich) galaxies. 

To date, metallicity measurements for high-redshift (z > 1) 
galaxies have been made via indirect methods from the ra- 
tios of strong nebular lines. In local galaxies, the preferred 
method for determining metallicity is via the "direct" Tg 
method which relies on the electron temperature sensitivity 
of the [O III] A4363 line. At high redshift, metallicity mea- 
surement via the direct method has not been possible because 
the [O III] A4363 line used to derive the electron temperature 
is very weak. The highest redshift t hat the [O III] A4363 line 
has been reported i s at z ^ 0.8 by iHoyos et al.l (l2005 h and 
iKakazu et all (l2007l) . 

Strong tensing by galaxy clusters offers an alternative tool 
to explore the metallicity of galaxies at high-redshift. One ad- 
vantage of tensing selected galaxies is that they are immune 
to color and luminosity biases. Previously, only five lensed 
galaxies had been observed with sufficient spectral coverage 
and S/N to yield metall icities ( Yee e t al. 1996; Pettini et al. 
[2000, 2002; Lemoine- Busserolle et al.ll2003l: iNesvadba et al.l 
12007: Stark et aL 2008) . Compared with Lyman break se- 
lected galaxies in the same redshift and luminosity range, 
the five lensed objects span a substantially broader metallic- 
ity range (12 + log(0/H)~7.1 — 9.1), indicating that current 

* Based on data collected at Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan. 



metallicity history studies at 2; > 1.5 (12 + log(0/H)~8.4 — 
8.8) may not be sampling the full metallicity range at this 
redshift. It is worth noting that the metallicity range quoted 
for the lensed galaxies is extremely uncertain, depending on 
various assumptions about reddening and metallicity calibra- 
tors. Another great advantage of gravitational tensing is its 
natural magnification of the source. The flux of the lensed 
background galaxy is commonly boosted by a factor of 1- 
3 mag or more around the critical lines of the cluster cen- 
ter. Optical spectroscopic surveys of high-z galaxies near 
cluster critical lines have already yielded a large catalog of 
lensed galaxies beyond z ^1 (e.g., Broadhurst et al. 2005; 
lRichard etal1l2008l: iFrve et al.l l20Q7l: iLimousin etaiTl2007h . 
Unfortunately, spectra from these surveys cannot be directly 
utilized in metallicity studies because the commonly used 
metallicity diagnostic emission lines have shifted into the near 
infrared (near-IR) at z > 1. 

With the revolutionary multi-object cryogenic near-IR 
spectrographs that have been recently installed on 8- 10m class 
telescopes, high efficiency near-IR spectroscopy is now avail- 
able. We carried out the first near-IR spectroscopic survey 
of lensed galaxies behind tensing clusters, aimed at measur- 
ing the metallicity of galaxies at 1.5 < z < 3. This Letter 
reports the rest-frame optical spectrum of a strongly lensed 
star-forming galaxy at z = 1.7 behind the Abell cluster 1689 
from our ongoing survey. The galaxy cluster Abell 1689 was 
chosen as our first target because it has the largest Einstein ra- 
dius of ~ 50^' dBroadhurst et al.l2005|) . and nearly ^ 100 spec - 
troscopically confirmed lensed images (ILimousin et al.l20Q7|) . 
The gal axy in question i s designated as ID 22.3 in the cat- 
alog of Broadhurst et al. ( 2005), we refer to it as Lens22.3 
hereafter. Throughout we use the standard lambda cold dark 
matter (ACDM) cosmology with Hq = 70 km s~^ Mpc~^, 
Qm = 0.3 and ft a = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

The data were obtained on 2008 March 26th, using the 
Multi-Obje ct InfraRed Ca mera and Spectrograph (MOIRCS) 
at Subaru (Ichikawa et al.h200 6). in the spectroscopic mode 
(Tokoku et al. 2006). We used the low-resolution grism zJ 
500 with a nominal resolution of 500 and a wavelength cov- 
erage of 0.9 — 1.78 fim. We designed two masks targeting 
the spectroscopically confirmed lensed images in the field of 
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Fig. 1 . — Observed Subaru MOIRCS spectrum of Lens22.3 at z= 1.7. The 
extracted one-dimensinoal spectrum is shown in black, with the identified 
emission lines labeled. The sky telluric absorption and OH emission lines are 
shown in red and blue, respectively. The two-dimensional spectrum for each 
of the detected lines is shown on top. 

Abell 1689. We obtained useable data for one mask. Dur- 
ing the observation for this mask, the seeing varied between 
(X'S and (y.'6, and the airmass ranged between 1.1 and 1.7. 
Lens22.3 was placed at one of the center slits of the mask, 
with a slit width of ~ t' and a slit length of ~ ^' . The 
telescope was nodded along the slit in a ABAB mode with 
a dithering length of 2'.'5. The exposure time for each individ- 
ual frame was 400-600 seconds, resulting in a total integration 
time of 6800 seconds. 

We reduced the data using the standard IRAF and IDL rou- 
tines. Sky OH lines were removed by subtracting the closest 
dithering pair A and B. A scaling factor between each pair 
was determined to optimize the sky subtraction. The sky sub- 
tracted A-B pair was flat-fielded, cleaned for comic rays, and 
then traced and rectified using a dome flat spectrum. The fi- 
nal two-dimensional spectrum was created by co-adding the 
reduced individual pairs. Wavelength calibration was carried 
out using the strong OH lines between 0.9 and 1.8 /im. A 
F5V standard star was observed at approximately the same 
airmass as the target. The Telluric absorption was obtained 
by dividing the star spectrum with a fitted blackbody curve. 
The resulting stellar spectrum was then divided into the ob- 
ject spectrum to correct for telluric absorption. 

3. ANALYSIS AND RESULTS 

3.1. Metallicity 

At redshift 1.7, all the rest- frame optical metallicity diag- 
nostic lines fall into the spectral coverage of MOIRCS zJ500 
grism. We have detected all the rest frame nebular emis- 
sion lines from [O ll] to Ha, including the usually very weak 
[Olll] A4363 line. Figure [T] shows our final two-dimensional 
and one-dimensional rest- frame optical spectrum of Lens22.3. 
Table [T] lists the measured line fluxes relative to H/3, the 
1 (J error, and S/N. Since the metallicity studies use only 
line ratios, we did not attempt to flux calibrate the spectrum, 
which is nontrivial in near-IR observations. The instrument 
response is essentially constant for the wavelength range from 
[O III] A4363 to [O III] A5007. The Ha line lies in the region 
where the response has dropped by ~ 35%. The response 
on the expected [Nil] line position has dropped by ~ 65%, 
which may be the reason that there is no detection of [N ll] 
within 3<j. After correcting for instrumental response, we ob- 



TABLE 1 

Observed Line Intensities and Ratios 



Ion Line 


^ohs 


Flux 


Flux (corrected) 


Icr 


S/N 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


[Nil] A6584 




<0.009 


<0.025 






Ha A6563 


1.705 


3.12 


5.03 


±0.4 


13 


[Oiii] A5007 


1.706 


6.45 


6.45 


±0.3 


15 


[Oiii] A4959 


1.705 


1.98 


1.98 


±0.3 


10 


H/3 A4861 


1.705 


1.0 


1.0 


±0.1 


5 


[Oiii] A4363 


1.696 


0.27 


0.27 


±0.1 


3 


[On] A3727 


1.708 


0.95 


1.11 


±0.3 


8 



Note. — (1) Detected emission lines and their rest frame wavelength. 
(2) Redshift calculated from the observed wavelength. (3) Measured flux, 
relative to H^^, detector response uncorrected. (4) Measured flux, relative to 
\{(3, detector response corrected. (5) 1 a error of the line flux from Gaussian 
fitting. (6) Signal-to-noise of the line flux. 

tain an upper limit for [N ll]/Ha. IVEetallicities were calculated 
using the following two methods. 

1. Strong line diagnostics. In the absence of direct mea- 
surements from electron temperature, the ratio of strong neb- 
ular lines is used to determine the oxygen abundance. The 
well-known discrepancies between the different line ratio in- 
dicators can be removed by conver ting to the same base cal - 
ibration following the methods of IKewley & EUisonl (|2QQ8|) . 
We use the indicator i?23 = ([Oil] A3727 + [Oiii] AA4959, 
5007)/H/3 to calculate metallicity. In our case, the upper and 
lower branch degeneracy of i?23 can be broken by the upper 
limit of [Nll]/Ha < 0.005. In addition, the detection of the 
weak [O ill] A4363 line almost certainly means the adoption 
of the lower branch. For t he actu al calculation, we use the 
lower branch calibration oflMcGaush (199 11) (M91 method ) 
in the analytical form given by Kobulni cky et al.l ('1999'). 
The resulting metalHcity i s 12 + log (0/H)= 8.0 in the M91 
method. Converted to the iKewlev & Dopltal (l2002h method 
(KD02), we have 12 + log(0/H)= 8.1 assuming zero extinc- 
tion, or 8.3 ass uming a high extin ction value of E(B-V) = 0.4 
(ICardelli et al.lll989l: lQsterbrocklll989h . The high extinction 
value agrees with the direct estimate from the in Ta- 
ble [TJ However, because the Ha line lies near the edge where 
the response has dropped by ~ 35%, the ratio of Ha/H/3 can 
be largely uncertain because it depends sensitively on the pre- 
cision of correction for the detector response. Therefore in- 
stead of using only the high extinction from Ha/H/3, we use 
two extreme values of extinction to quantify the range of un- 
certainty caused by extinction. Note that the logarithm in the 
definition of 12 + log(0/H) means that 10% of change in the 
line ratio causes only ~ 0.04 dex change in metallicity. The 
uncertainty for our derived metallicity is therefore < 0.1 dex 
prior to extinction correction. 

2. Direct metallicity from electron temperature Tg. The 
auroral line [O III] A4363 is often very weak, even in low 
metallicity environments. It is usually not observed in high 
metallicity galaxies without very sensitive, high S/N spectra 
(e.g., Garnett et al.l l2004D . The 3 cr detection of [O III] A4363 
in Lens22.3 strongly suggests a low metallicity environment. 
The ratio of [O ill] A5007, 4959 and [O ill] A4363 allows a di- 
rect measurement of the oxygen abundance via electron tem- 
perature Tg. We use the nebular package in IRAF to de- 
rive the electron temperature Te = 2.3557 x 10^ K, assum- 
ing an electron density ng of 100 cm~^ (Tg is relatively in- 
sensitive to ne, e.g., for Ue = 1000 cm-^, Tg = 2.3508 x 10^ 
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K). Following the procedure of llzotov et alJ (l2QQ6h . we ob- 
tain 12 + log(0/H) = 7.5 ±0.1 assuming zero extinction; or 
12 + log(0/H) = 7.3 ± 0. 1 assuming a high extinction value of 
E(B-V) = 0.4. 

The 0.6 dex discrepancy between the KD02 and Tg method 
is not surprising since the metallicity measured from the di- 
rect Tg method is systematically lower th an other methods ; 
such offsets can be as large as 0. 7 dex (iLiang et al.l l2006t 
lYin etal.ll2007l: iNagao et al.lf2006h . Regardless of the off- 
set between the two methods, both metallicity estimates in- 
dicate that Lens22.3 is low in oxygen abundance. The ev- 
idence that Lens22.3 is a star-forming galaxy instead of a 
narrow line AGN comes from its position on the traditional 
Baldw in — Phillips — Terlevich (BPT) diagram (Bal dwin et al.l 
Il98l|) . Using the upper limit of [Nll]/Ha and the ratio of 
[Olll]/H/3, Lens22.3 is located on the leftmost region of the 
star-forming branch of the BPT [N ll] diag ram, e xcluding its 
possibility of being an AGN ( Kewlev et al.ll2006l) . 

3.2. Photometry 

We searched the Hubble Space Telescope (HST) archive 
and obtained the photometry for Lens22.3. The observed AB 
magnitudes are 23.2 ± 0.1, 23.6 ± 0.2, 23.3 ± 0.1, 23.4 ± 
0.1 and 23.1 ± 0.1 in the ACS gpAibW^ ^f625T^, 'ipiibw^ 
ZF850LP, and NICMOS J mow bands, respectively. After 
correcting for the lensing magnification of 15.5 ± 0.3 (pro- 
vided by J. Richard 2009, private communication), the abso- 
lute B band magnitude is = -18.3 ± 0.1 (Vega System). 
A simple spectral energy distributi on (SEP) fitting using 
the kcorrect v4_l_4 code of Blanton & RoweisI (l 2007h 
yields a stellar mass of 1.6 x 10^ Mq. Since the kcorrect 
may overestimate the stel lar masses of high-redshift galaxies 
(|Blanton & RoweisI |2007|) . we recalculated the stellar mass 
using the advanced version of the code LE PHARE (pho- 
tometric redshift and simulation package; Ilbert et al." "2009b 
based on population synthesis models of Bruzual & C harloi 
(l2003h . The resulting stellar mass is 4.4±1.2x 10^ M©. 

3.3. Morphology and Kinematics 

Thanks to the excellent seeing (~ 0'.'4) and the powerful 
magnification of gravitational lensing, the target is spatially 
resolved. In Figure [21 we show the HST ¥115^ image in 
the vicinity of Lens22.3. The image has been registered onto 
the MOIRCS chip, which has a pixel scale of 0^.^117 pixel" ^. 
The MOIRCS slit layout is shown, where X is the dispersion 
direction and Y is the spatial direction which is also along 
the elongated side of the galaxy. On the lower right corner 
is a zoomed image of Lens22.3, with various intensity con- 
tours overlaid. It has a clear double-nuclei-like morphology. 
The nuclei are separated by ~ 0'.'4, corresponding to a phys- 
ical scale of ~3 kpc at the redshift of z = 1.7. The intrinsic 
separation will be smaller than this value because the image 
has been elongated by lensing effects. We divide the two- 
dimensional spectrum into subsections A, B, C, D along the 
spatial direction Y, using a 3 pixel bin that roughly matches 
the average seeing of the observation. We extract the spectrum 
for apertures A, B, C, D using the strongest line [O III] A5007 
and Ha. The velocity profile relative to the center of the sub- 
sections are given in Figure [3] For comparison, the velocity 
profile for the sky OH line nearest to Ha extracted from the 
same apertures is also shown. There seems to be a systematic 
rise in the velocity from A to D in Ha, but not in [O III] A5007. 
The misalignment of the line profile in the spatial direction 
and the double-nuclei like morphology may indicate that the 
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Fig. 2. — //^r rescaled F775W image for the field of Lens22.3 The image 
has been rescaled and registered onto the MOIRCS chip. The slit on Lens22.3 
has a width of V and a length of Slit X (dispersion) and Y (spatial) 
direction are marked. A zoomed view of Lens22.3 with intensity contours 
(0.3, 0.2, 0.1, 0.04 of the peak intensity) is shown on the lower left corner. 




Velocriy (Km/s) 

Fig. 3. — Top two panels: velocity profile of Ho; and [Oiii] from spectra 
extracted within the spatial apertures A, B, C, D (bottom right panel). Bottom 
left panel: velocity profile for the sky line nearest to Ho; using the same 
apertures. Bottom right panel: two-dimensional spectrum of Ha, spatially 
cut into sub-apertures A, B, C, D. The velocity zero point is chosen to be the 
averaged center of the cuts. 

system is rotating or merging. Future follow up with integral 
field spectroscopy with adaptive optics is highly desirable to 
reveal the kinematics of the system. 

3.4. Mass -Metallicity Relation 

In Figure m we sh ow the mass/luminosi ty — metallicity re- 
lation f or the local (iTremonti et al.l [2004|) and high-redshift 
(lErb et al. 2006) star-forming galaxies. All metallicities in 
Figured] have been rescaled t o the s ame KD02 method using 
th e calibration of Kewlev & EUisonl ([2008). The local sample 
of ITremonti et al.r(l2004l) is based on - 53,000 Sloa n Digital 
Sky S urvey (SDSS) star-forming galaxies and the Erb et al.l 
(l2006h sample is based on composite spectra of 87 rest-frame 
UV- selected star-forming galaxies at z > 2. Lens22.3 is con- 
sistent with the mass/luminosity — metallicity relation of the 
high z sample. We evaluate the effect of extinction assum- 
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Fi g. 4. — Left: mass — metallicity relation. Dash line: local relation from the SDSS sample ("Tremonti et al."2004'). Filled black circle: relation for galaxies z 
> 2 dErb et a l. 2006). Filled red square: Lens22.3, without extinction correction. Empty red square: Lens22.3, after extinction correction. Right: luminosity — 
metallicity relation. Dash line: local relation from SDSS sample ( Tremonti et al- 2004). Filled black circle: relation for galaxies z > 2 (Erb et al. 2006). Filled 
red square: Lens22.3, without extinction correction. Empty red square: Lens22.3, after extinction correction. All metallicities have been recalculated to the same 
base of KD02 ( Kewley & EUison 2008). 

ing a high extinction value of E(B-V) = 0.4. Extinction causes 
metalHcity to rise by ~ 0.1 dex, but it does not significantly 
change the position of Lens22.3 on the mass — metallicity di- 
agram. Lens22.3 is currently the lowest mass system with di- 
rect measurement on the mass — metallicity relation diagram; 
Lens22.3 extends the luminosity — metallicity relation to the 
fainter end by more than an order of magnitude, demonstrat- 
ing the capability of using gravitationally lensed samples to 
probe intrinsically fainter and less massive systems than are 
found in existing samples. 



luminosity-metallicity relation by more than an order of mag- 
nitude. Our on-going near-IR spectroscopic survey of lens- 
ing clusters aims to obtain metallicities for more objects like 
Lens22.3 in the redshift range 1.5 < z < 3. 



4. SUMMARY 

We reported the discovery of a low metallicity star-forming 
galaxy at z = 1.7. We detect the temperature sensitive au- 
roral line [O III] A4363 which allows us to derive metallic- 
ity using the direct Tg method for the first time at high red- 
shift. The low metallicity and low luminosity /mass of the 
z = 1.7 star-forming galaxy place it at the lowest end of the 
mass/luminosity — metallicity relation, extending the z > 2 
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